Reconstruction of superoperators from incomplete measurements 
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We present strategies how to reconstruct (estimate) properties of a quantum channel described 
by the map £ based on incomplete measurements. In a particular case of a qubit channel a complete 
reconstruction of the map £ can be performed via complete tomography of four output states £[qj\ 
that originate from a set of four linearly independent "test" states Qj {j = 1, 2, 3, 4) at the input of the 
channel. We study the situation when less than four linearly independent states are transmitted via 
the channel and measured at the output. We present strategies how to reconstruct the channel when 
just one, two or three states are transmitted via the channel. In particular, we show that if just one 
state is transmitted via the channel then the best reconstruction can be achieved when this state is 
a total mixture described by the density operator g = 1/2. To improve the reconstruction procedure 
one has to send via the channel more states. The best strategy is to complement the total mixture 
with pure states that are mutually orthogonal in the sense of the Bloch-sphere representation. We 
show that unitary transformations (channels) can be uniquely reconstructed (determined) based on 
the information of how three properly chosen input states are transformed under the action of the 
channel. 
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I. INTRODUCTION 



Our observation of external world is always incomplete 
and the main task of a physicist is to find the best possible 
description of our state of knowledge about the physical 
situation. A typical example of this is a description of the 
action of an unknown quantum channel £. Specifically, it 
is well known that a complete correlation measurement 
between a properly chosen set of linearly independent in- 
put states pj and corresponding outputs £[pj] reveal the 
complete information about the character of the map £ 
(see e.g. Ref. [2-4]). The question is, how to describe the 
map £ if only results of incomplete correlation measure- 
ments are known? In this paper we will address some 
aspects of this problem. 

Let us formulate our problem in more detail: In quan- 
tum theory (under reasonable conditions) the evolution 
map is associated with a trace-preserving completely pos- 
itive linear map £ [5,6] acting on a set of quantum states 
S{H). The set of states is a subset (not subspace) of the 
real linear space of all hermitian operators with dimen- 
sion (P, where d = dimTi. To define a linear map £ we 
have to specify its action on an operator basis, i.e. we 



true physical situation ... is 

unknown, and always remains so, 

because of incomplete information. 

E. T. Jaynes [1] 



have to specify d^ x d^ parameters [2]. In addition, the 
condition that the map £ is trace-preserving reduces this 
number to d^ — d^ . From here it follows that in order 
to determine a transformation £ experimentally we have 
to choose d^ linearly independent states pj . These states 
are sent via the channel £. At the output of the channel 
we perform a complete state reconstruction of the states 
£[pj\. Given the fact that each of the d^ states is de- 
termined by d^ — 1 real numbers we clearly see that the 
map £ can be completely determined via this so-called 
quantum process tomography. 

However, also another scenario is possible [3,4,7,8]: 
One can use a single maximally entangled state to probe 
the action of the channel. In particular, let us consider a 
maximally entangled state ]V'+) = -j^^j \j) ® \j) from 
the Hilbert space Ti^Ti. One of the subsystem is trans- 
formed by the channel, while the other (reference) system 
evolves freely without any disturbance. In this case the 
final state of the composite system reads 



n = £(^I[P+] = ^^f[ejfc]' 



' ejk ■ 



(1.1) 



jfe 



where P+ — ]-0_|-) (?/'+] = ^ J^jk ^jk ^ Cjk and operators 



*We dedicate this paper to Asher Peres on the occasion of his 70th birthday 



Cjk — \j){k\ form an operator basis of the operator space. 
Via the reconstruction of the output state f2 we obtain 
the complete information about the map £. This is due to 
the fact that the input-output correlation between states 
P+ and Q provides us with the knowledge how the oper- 
ator basis is transformed. 

In what follows wc will consider the first reconstruction 
scenario (the one that does not involve entangled states) 
but we study the situation when the number of linearly 
independent input states is less than cP. In this situation 
the perfect process tomography cannot in general be per- 
formed. To simplify our analysis we will assume that the 
preparation of test states pj as well as the determination 
of the output state can be performed perfectly (e.g., via 
quantum-state tomography) . 

The present article is organized as follows: In Sec. II 
we briefly review basic formalism of completely positive 
maps that act on a qubit. For the completeness of our 
discussion we will also show how the CP map can be 
completely reconstructed. In Sec. Ill we will present our 
strategy how to perform reconstruction of quantum pro- 
cesses from incomplete measurements and we will analyze 
number of examples. We conclude our paper by several 
comments that can be found in Sec. IV. 



II. PROPERTIES OF QUBIT QUANTUM 
OPERATIONS 

In what follows we will focus our attention on a recon- 
struction of single-qubit transformations - i.e., CP maps 
on qubit channels. The state space S{H) is a convex sub- 
set (not subspace) of a four-dimensional linear space of 
hermitian operators. This means that any state can be 
expressed as a four-dimensional real vector. Due to lin- 
earity of quantum mechanics, the state transformation £ 
can be represented by a 4x4 matrix. If we use the basis 
of a matrices (ctq = /, cti = ax, 12 = (JyiCr^ = az), then 
an arbitrary state of a qubit can be expressed as 



'=2^' 



^) = 



1 



f= {x,y,z), 



(2.1) 



with x^ + j/^ + z^ < 1. From here it follows that 
the state space of a qubit can be represented (visual- 
ized) as a three-dimensional (Bloch) sphere. The center 
r = (0,0,0) of the Bloch sphere corresponds to an equally 
weighted superposition of all qubits states, i.e. it repre- 
sent a total mixture. Let us note that the Bloch sphere 
(state space) form a subset of a three-dimensional space 
of hermitian matrices with unit trace. In what follows 
we will often omit the factor 1/2 in the vector represen- 
tation. 

In this basis any evolution £ takes the affine form 



£ 



1 
t E 



(2.2) 



where E is a, real 3x3 matrix. In particular, 

r ^t r = Er + t . 

The vector t corresponds to a translation of the center 
of the Bloch sphere. If f = (J, then the total mixture is 
preserved and the map £ is called to be unital, or equiv- 
alently bistochastic [6] . For a unitary transformation the 
matrix E corresponds to a three-dimensional orthogonal 
rotation, i.e. E~^ = E^ . The Bloch sphere picture en- 
ables us to illustrate quantum operations as deformations 
of the unit sphere [6,9,10]. 

The matrix form (2.2) of £ guaranties the condition 
of the trace preservation (first row) . As we have already 
mentioned the matrix £ is specified by 2^ — 2^ = 12 real 
parameters, but because of the constraint of the complete 
positivity the range of these parameters is restricted. The 
transformation £ is completely positive if and only if it 
can be written in the Kraus form [6] 



m 






AkgAi , with 



E4^^ 



(2.3) 



or equivalently [11] if the operator 51 = £ ® X[P+] is pos- 
itive and if Trifl = 1/2, where Tri corresponds to a par- 
tial trace over the first system. For completely positive 
maps the following relations hold [6]: 



1. Contractivity of distance 



D{g,0>D{£[g],£[^]). 



(2.4) 



where D{A,B) = Ty\A - B\. For qubits this con- 
dition reads 



D{g,0 = \r-^, 

where g ^^ r and ^ <-> s . 
2. Fidelity monotone 

F{g,0<F{£[g],£[^]) 



(2.5) 



(2.6) 



where F{A, B) = TrVvABVS is the fidehty. 



For our purposes it is useful to present explicitly the 
distance and the fidelity between any state g and a total 
mixture ^I ^^ rn = (0,0,0), i.e. the center of the Bloch 
sphere. This distance equals 



Dig,^!) 
1 



1 



1 



F{g, -I) - ^TrV^ = 3 (v/T+M + ^1^=^) • 

Both of these equations are monotone in the state g. The 
fidelity decreases as the |r| increases and the trace dis- 
tance increases together with \f]. For non- unital maps 



(f [^/] 7^ |/) there are states for which D{q^ i/) < 
D{S[q], i/), or F{q, \I) > F{£[q], i/). Therefore if we 
find such states then we can conclude the map £ must 
be non-unital. 

Any state transformation £ can be (uniquelly) ex- 
pressed as a composition of maps [9] 



£[g] ^UVV[q\ = UV[VqV'^]U^ 



(2.7) 



where C/, V are unitary operations and I? is a map of the 
form (2.2) with a diagonal matrix D — diagjAi, A2, A3}. 
The A's are singular values of the matrix E associated 
with the original map £, i.e. the squares of the eigenval- 
ues of the matrix EE^ . We should note that also the vec- 
tor t is changed according to the rule t -^ Rjjt where Ru 
is a rotation corresponding to the unitary transformation 
U . The existence of such representation follows from the 
polar decomposition theorem. According to this theorem 
matrices E can be expressed in the form E = RjjDRv 
where Rjj^ Rv are orthogonal rotations. This means that 
the conditions due to the complete positivity of the map 
£ are reduced to the requirement of the complete posi- 
tivity of the map T). Unfortunately, even for a qubit the 
conditions of complete positivity of T) is quite complex. 
The necessary conditions are [12,9] 



{\i±\2f<{l±\^f-tl 



(2.8) 



If ti — t2 = 0, then these conditions are also sufficient 
ones. Consequently for unital maps {t = 0) we have the 
necessary and sufficient conditions in a simple form 



(Ai±A2)'<(l±A3 



(2.9) 



The set of all completely positive linear trace- 
preserving maps (quantum channels) is convex. That 
is, we can define "pure" quantum channels as extremal 
points of this set. Obviously, unitary transformations are 
extremal. However, there exist also non-unitary extremal 
quantum operations [9]. For instance, a contraction to 
single pure state, i.e. g -^ |V')(V'I for a-U Qj is extremal. 

In what follows we will often use the following basis of 
the state space 



In what follows we will use the operator basis {Sj} which 
is obtained by a unitary transformation of the original 
cr-basis. This basis can be illustrated as a set of orthog- 
onal axes in the Bloch-sphere picture. The above basis 
is determined by the choice of cartesian coordinates in a 
three-dimensional space. 

With the help of the four states gm, fe, gy, gz we can 
perform the complete reconstruction of the map £. The 
reconstruction of £ requires the knowledge of the trans- 
formation of these states, i.e. f [pm],^ [Pa;],^ [g'l/],^ [giz]- 
These output states determines columns of the matrix 
£. In particular, using the correspondence £[gx] ^^ x' 
(similarly for other basis operators) we can write t — m! 
and the matrix E reads E = {x' — t^if — t,z' — i), so the 
transformation £ can be expressed as 



£^ 



1 















m X — m y — m z — m 



(2.11) 



In general we can use any four linearly independent states 
for reconstruction purposes. Nevertheless, it should be 
noted that it is always possible to transform the recon- 
struction procedure into the basis discussed above. 

As we have already mentioned, instead of using four 
linearly independent states of a qubit, one can use a sin- 
gle entangled stateP+ of two qubits to perform a perfect 
reconstruction of the channel £. Specifically, from the 
output state ^ = £ ® 1[P+\ we can determine the map 
£. 

In the (J basis the state P+ reads 

P+ = -{1 ® 1 +ax®(Jx- (jy®ay + (j^® a^) . (2.12) 
The output state g can be uniquely written in the form 



n-- 

where 



-{^o®l 



no 



■ ^x (Xi 

-£[1]- 
-£[a-x] 



i i-ji 



n. 



1 



i U'jj Cy ^U ij 



2^^ 



Qrn 



gy = 



; Qz = 



These states are eigenstates of the operators Sx,Sy, 



Sz 



associated with the value -1-1. The state m describes 
the total mixture. The operators Sj are the generalized 
CT- matrices with respect to a qubit state lip). They are 
defined by the following relations 



5, = -z(|^)(V'i.|-|V±)(V'l), 
Sz^\i^){i^\-\i^±){i^±\. 



(2.10) 



^z =£['^z] = 7fz-d . 

The reconstruction of the state O gives us the operators 
rife and the map £ is completely determined by the above 
equations. In particular, the vectors e^ (with cq = i) cor- 
respond to columns of the matrix £, i.e. 



£ 



10 

X G-x Gy €z 



(2.13) 



Before we conclude this section we note that in general, 
we do not need maximally entangled states to perform 
the complete reconstruction of the map £ (for details see 



Ref. [3,4]. In fact, any state that can be written in the 
form 

g= -{1 ® 1 + xcFx ® (yx+ yay ® ay + za^ ® ctz) (2-14) 

can be used for the reconstruction, because the output 
state reads 

fi = - {£[1] 'S) 1 + x£[(73:] ig) ax + y£[fTy] <8) CTy + z£[az] <S) az 

and the one-to-one correspondence between £ and fl is 
obvious. In particular, every pure state can be written 
in the Schmidt form \ip) = ajoO) + 6|11), where {|0), |1)} 
are suitable chosen bases of subsystems. If we define the 
operators S^ = |0)(1| + |1)(0|, Sy = -i|0)(l| + i|l)(0|, 
^, = |0)(0|-|l)(l|,thcn 

n= -{1 <» 1 + 2abSx ®Sx~ 2abSy <» Sy + Sz<» &) . 

Consequently, whenever the pure state is correlated (en- 
tangled), it can be used to perform the complete recon- 
struction of the CP map. 



III. STRATEGIES FOR INCOMPLETE 
RECONSTRUCTION 

Any reliable reconstruction scheme should be as con- 
servative as possible, i.e. we cannot prefer any type of 
transformations that is not warranted by the data [1]. 

Having this principle in mind let us try to guess the 
expression for the map (i.e. the transformation) when 
no measurement is performed. In the absence of mea- 
surement data arbitrary test state at the output of the 
channel can be estimated as a total mixture described by 
the density operator 1/2 (for more details sec [1]) Conse- 
quently, the estimated map describing the action of the 
channel describes a contraction of the whole Bloch sphere 
into a single point - the total mixture, i.e. Eq : g i-^ ■^I . 

The basic idea of our strategy is simple. If wc have no 
prior information about the action of the channel and the 
data do not contain any information about the transfor- 
mation of a state g, then we assume that any input state 
is transformed by the channel into the total mixture, i.e. 
£[q] — 2-^' However, we have to clarify one point: If we 
know a priori that the total mixture is affected by the 
channel £, i.e. i?[5/] 7^ 5/, then we cannot assume that 
the unknown input g is mapped into the total mixture. 
We have to assume that it is mapped as £[g] = i£[^/]. 
Of course, our guess must satisfy constraints imposed by 
the complete positivity. Let us now present our strategy: 



2. Step 2. If we have no information about a "shift" of 
the total mixture, then we must check whether the 
map £ is unital providing that all unknown states 
are mapped into ^/. In other words, check the 
complete positivity of the map £ when all free pa- 
rameters are equal to zero. If yes, our guess is com- 
pleted. 

3. Step 3. Choose the free parameters such that £ is 
I completely positive, D{£[^I],^I) is minimal and 

D{£[gj\,£^i^[-^I]) is also minimal for all undeter- 
mined states g. 

In what follows we will analyze reconstruction of qubit 
channels. We will study three relevant cases: When we 
know how a single state is transformed, when we know 
how two (linearly independent) states are transformed 
and finally, when we known how three states are trans- 
formed. 



A. No input state 

Our guess is very simple in this case. Our strategy im- 
plies that all inputs are mapped onto the total mixture. 
So the whole Bloch sphere is contracted into its center 
and the map reads 



£0 = 



10 









(3.1) 



B. Single input state 

In this case we have to consider three different cases. 
Either the input state is the total mixture, or it is a pure 
state, or an arbitrary mixed state. 

Let us assume that our information is the following: 



w 



■ (?) , that is we know that the input state 



that has been prepared in a total mixture has been trans- 
formed into the state ^{I + 7^ ■ a). In this case our guess 
of the transformation of other basis operators is given by 



£[qx] =f [^y] ^£[qz] 



.(!., 



1 



(I 



Consequently, the reconstructed quantum channel £ is 
the contraction to the point £[\l] which is obviously a 
completely positive map. 

Consider the case, when we know how a pure state \ip) 
is transformed, i.e. 



ma^i 



1 



1 



1 



{I + Sz)^-(I + r' -3) = -{I + r'S',) 



1. Step 1. Check whether the available data contain 
any information about the transformation of the 
total mixture, i.e. whether we know f [j-^]- 



where r' = \f'\. We can assume that the map £ is unital, 
because D{£[P^], \l) = r' < I = D{\l,P^). The guess 
is that £\gx\ = £\Qy\ ~ ^[2^] ^ 2^ transforms the whole 



Bloch sphere into a section of a line containing the to- 
tal mixture as its center. Note that Qx^ Qy are expressed 
with respect to the operator basis with Sx,Sy,Sz, i.e. 
\4'){'4'\ = Pi> = Qz- In fact, such transformation is com- 
pletely positive. The matrix £ takes the form 



£ = 



10 
7=' 



(3.2) 



The singular values of the matrix E are equal to Ai — 
A2 = and A3 = r'. Consequently, the sufficient crite- 
rion for the complete positivity of unital maps (2.9) is 
fulfilled. 

Finally wc shall discuss the case when we know that the 
map under consideration performs the following transfor- 
mation of a single input state 



1 



g=-{I + rS,) 



1 



(I + r'S'J 



(3.3) 



In this case we must investigate whether the map £ can 
be still unital, or whether the total mixture is shifted 
from the center. One can easily see that 

so the unitality depends on the difference A = r — r' . 
li A > {g IS shifted closer to the center), then the 
transformation might still be unital, but for A < the 
quantum channel is inevitable nonunital. 

The first guess can be, that we should move the total 
mixture as little as possible so that D{£[g],£[^I]) — r, 



Sill] 



{I + {r'-r)Si) 



(3.4) 



In this case D{g',£[^I]) — \(r' — r) — r'\ = r. Unfortu- 
nately, the assumption (3.4) about £[^I] is not compati- 
ble with the positivity condition. Note that a linear (not 
convex) combination A = Xg + {1 — X)^I = i(/ + XrSz) 
for A = 1/r corresponds to a pure state gz = ^{I + S^)- 
Applying the knowledge about the state transformation 
£ we obtain 



£[qz] ^ -g' 

r 



r/ ^2 ^ 2 ^ 



(1-A)5^) 



with A = r — r' < 0, i.e. the operator £[gz] is not posi- 
tive, an therefore the map i?[^/] cannot be defined in this 
way. So the shift must be such that also the operator 
A = ^£[g\ + {1- j:)£{^I] is positive, i.e. A = ^{I + nS'z) 
with K < 1. Comparing two expressions for A we obtain 
that 



4'^-b' 



1 -r 



-Si 



The distance i:)(£[i/],i/) = |(r'-Kr)/(l-r)| is minimal 
for K = 1 (under the condition that k < 1). Consequently 
the guessed state transformation is 



Qxt gyi n^ 



4'14 



1-r 



-s: 



We obtain that the pure state gz is mapped always into 
a pure state g'^. It is easy to check that the map 



/ 1 





V ^ 1 

\ 1 — r 



\ 



r —r 
1-r 



(3.5) 



is indeed completely positive. Note that this is in fact 
the expression of the matrix T) associated with the map 
£, i.e. it is not an expression of the matrix in a fixed 
basis, but rather a matrix that relates the two bases 
{SxTSy^Sz} and {S'x,Sy,S'z}. Our data contain infor- 
mation only about the operators Sz and S'z- The op- 
erators Sx,Sy,S'x,S'y are arbitrary, such that TrSjSk — 
TiS'jS'f. ~ 6jk, i.e. in the Bloch-sphere picture they define 
three mutually orthogonal axes (a cartesian coordinate 
system). 

The case when r' < r (A > 0) is easier. Then we can 
make a following guess 



g^ g 



Qx} gyj „ ^ 

and consequently the map 



.,1.1 



1/ 

2 



10 





r'/r. 



(3.6) 



is completely positive. Note that the matrix in Eq.(3.5) 
for r' = r takes the form (3.6). 

In general, when we do not know how the total mix- 
ture is mapped, our guess can be written in a compact 
form 



/ 



£^ 



1 















\ 





Vmax{r:^,0} min{i^, 1} - max{r^,0} / 

Note: Adaptable basis 

Our reconstruction results with the expression of the 
quantum operation that is not written in some fixed ba- 
sis. The matrix £ acts on the input state written in one 
basis {Sz,Sy,Sz}, but it produces output state written 
in the basis {Sx,Sy,S'z}- Our data determine only op- 
erators Sz, S'^. The choice of other basis operators must 
be made so that they are mutually orthogonal. How- 
ever, the choice of S'x, S'y is irrelevant, because our map 
transforms the whole Bloch sphere only into the subspace 
spanned by the identity / and the operator S[,, i.e. into 
the one-dimensional space (line) in the Bloch sphere pic- 
ture. The transformation £ can be expressed in the fixed 



basis just by applying a suitable rotation Rjj, i.e. the 
unitary transformation either before, or after the map £ 
is applied. 



1. Examples 

1. Identity. 

Let us assume that we find from an experiment that 
the single input state g is not affected by the channel. 
That is, the channel acts on the state g as an identity 
X and Q ^^ g, i.e. r' = r. In this case our guess of the 
channel is 



10 





1 



and under the action of this map the Bloch sphere trans- 
forms into a line. 
2. Unitary transformation. 

The case of unitary transformations is basicly the same 
as the case of identity, i.e. our guess is represented by the 
same matrix. However, we have also information about 
the transformation of one operator basis element, namely 



1 



Q^7;il + S.) 



where SL = f' ■ S with 



g' = UgU' = ^{I + S',), 



1. From here it follows that 



the quantum map £ can be expressed in the fixed basis 



10 

a;' 

y' 

z' 



(3.8) 



,z'). 



where we used a notation r' ~ (x', y' 
3. Contraction to a pure state. 

Let us assume that from an experiment we know that 
a test state is transformed by the action of the channel 
into the state P^: 



Q 



P, 



^■ 



(3.9) 



Let us consider that g is not pure. Then our guess al- 
ways concludes that the map is non-unital. Moreover, 
the matrix £ reads 



10 





10 



(3.10) 



so that all input states are mapped into the state P^. 
It means that such contraction is perfectly specified 
(guessed) just based on a single state transformation. 



Moreover, we have no doubts that the map £ has to be 
a contraction. 

The situation is different, if g is a pure state. Then the 
estimation (guess) of the map is expressed by Eq. (3.7). 
That is, the Bloch sphere is mapped into a line instead 
into the single point. 
4. Specific example. 

In this example we will apply our strategy to a partic- 
ular channel which in the cr-basis reads 



10 

0.5 0.2 -0.1 0.1 

0.2 -0.3 

0.3 0.3 



(3.11) 



(3.7) The action of this map is depicted in Fig.l. 




FIG. 1. The transformation of the Bloch sphere under the 
action of the map £ given by Eq.(3.11). The grid on the Bloch 
sphere and the ellipsoid helps us to see the rotation of the 
original state space. We also see that the map is non-unital. 

Let us assume that from the experiment we know how 
a single state g = 5 (-^ + 0.6fTa;) is transformed under the 
action of the transformation (3.11), i.e. 



1 



0.6a,) 



1 



(/ + 0.62cr, + O.Uay 




FIG. 2. The transformation of the Bloch sphere under the 
action of the reconstructed map £0, i.e. the contraction of 
the sphere into a single point representing the total mixture. 
Compare with the image of the real map £ in Fig. 1. 



It is easy to verify that given this transformation the 
total mixture has to be shifted, i.e. the (unknown) chan- 
nel is non-unital. Simply the inequality 



= v/O.622 + 0.122 > ^/oj2 



(3.12) 



implies the non-unitality of the transformation. The rea- 
son is that the state g' lies further from the total mixture 
than the state q. 

In the adaptable basis (see the discussion above) our 
guess has the form (3.5) 



^ad 



/I 

0.079 0.921 



V 



(3.13) 



This matrix transforms between two bases: {Sj — CTj} 
and {S', = ^{Oma^ + Q.l2ay),S'y = i(0.12a, - 
0.62o-y),S'^ — a^}, i.e. the output state £°''[gi] is ex- 
pressed with respect to the prime basis {S'^, Sy,S'^}. In 
the cr-basis the reconstruction reads 



£i 



1 

0.0776 0.904 

0.015 0.175 





(3.14) 




FIG. 3. The transformation of the Bloch sphere under the 
action of the reconstructed map £i given by Eq. (3.14). The 
sphere is transformed into a line. Compare with the image of 
the real map in Fig. 1. 



C. Two input states 

Let us assume that from an experiment we know how 
two input states are transformed under the action of an 
unknown map £. What would be our most reliable re- 
construction of the map? 

By knowing how two linearly independent input states 
are transformed by an unknown map £ we fix 6 parame- 
ters of this map. The other 6 parameters have to be spec- 
ified (guessed). We will follow our strategy described in 
the previous section and we will optimize our guess over 
six unknown parameters. 



The state transformation that are supposed to be 
known are: 



gi^ gi. 



and 



02^ g2 



(3.15) 



Thus, we have the information how a two-dimensional 
subspace oi S{Ti.) is transformed, or equivalently the one- 
dimensional subspace (line) in the Bloch-sphere repre- 
sentation. We can define two different operators (linear 
combinations of gi and 52) with more suitable properties 
for the reconstruction. These new operators are not nec- 
essarily positive but they have a unit trace. In terms of 
these new operators we obtain two new transformations 



Ci 






?i = 2'-^"''"'^^' 



(3.16) 



bS„ 



^2 = 7:iI + Psm9S'^+(3cosdS'), (3.17) 



that contain the same information about the map £ as 
the original ones. Transformation given by Eq. (3.16) 
specifies how a pure state fi is transformed, while 
Eq.(3.17) indicates how the "orthogonal" (see Fig. IIIC) 
operator ^2 is transformed. 

Let us consider firstly, that the line connecting gi and 
P2 contains a total mixture (^0 = 9-^)' ^-^^ ^^ know how 
the total mixture is transformed plus we have a knowl- 
edge about the transformation of one pure state f 1 . This 
is the same situation that had been considered in the 
previous subsection (now we are certain how the total 
mixture is transformed). Following our strategy we put 



^0 



1 



.,1.1 



1 



{I + Psm9S'^ + (3coseS')] 



2 



2' 
^2 ^ 1;{I + Sy) ^ £[ll] 



1 



e3 = ^(/ + 5,)^f[i/] = i(/ + /3sin<+/3cos05;), 
thus the map £ reads 



1 

/3sin 61 a-(3sin0 

/3cos6' -Pcos9 





(3.18) 



So under the action of this map the Bloch sphere trans- 
forms into a line. This map is completely positive only if 
a2-h4/32-4a/3sin6l< 1. 

Here we have to deal with the self-consistency of 
the available data, i.e. the transformations given by 
Eq. (3.15). Specifically, not all transformation gi —f g\^ 
g2 ^^ g'2 can be simultaneously associated with com- 
pletely positive maps. The sufficient and necessary con- 
ditions were derived in Ref. [13]. If for all positive t the 
inequality 



D{gi,tg2) > D{g[,tg'^) 



(3.19) 



holds, then there exists a completely positive map £ such 
that g'j = £[qj] (j = 1,2). Consequently this condition 
must be satisfied by our data. One can also find simpler 
necessary conditions that must hold. For instance the 
contractivity of the distance (if t = 1) of the map £ has 
to be fulfilled. 

In the case under consideration [see Eq. (3.18)] the con- 
dition D{^I, Qx) > D{£[\I],£[qJ) implies that q;^+/3^ — 
2a(3sm6 < 1. Obviously, also o? <\ and 0^ <\ must 
hold to ensure that f [^/] and £\q^ are density operators. 

If from Eqs. (3.15) we cannot infer the information 
about the transformation of the total mixture, then the 
estimated map takes the form [see Eqs. (3. 16-3. 17)] 



£^ 



10 

X a — X (/3sin0 — a;)/6 m 

y —y {pcos9 — y)/b n 

z ~z ~z/b k 



(3.20) 



with six free parameters x,y, z,m,n,k. Our strat- 
egy is to minimize the shift of the total mixture 
(i.e. y^x'^ + j/^ + y"^) and then minimize the distance 
D{£[^I],£[qz]) — \Jrn? + n"^ + k'^ under the condition 
that the matrix £ represents a completely positive map. 
In principle this task can be performed numerically, 
once we have particular values of b, a, (3, 6. Note that 
Qz = ^{I + S'z), where S'^ are determined uniquely once 
S'^,S' are given. Our data contain information about 
Sx, Sy, S'^, S' . The existence of £ is not guaranteed for 
all values oi a,l3,6. 

Let us first study the case, when it is possible to set 
f [i/] = /, e.g. the complete mixture is not affected by 
the map. In this case the transformation £ reads 



£ = 



/I 

a f3sm9/b m 

/3cos9/b n 

\0 k 



Now we have to determine the range of the parameters 
involved in this expression for which the transformation 
£ is completely positive. The contractivity condition of 
£ implies that our data must satisfy the relations 



a < 1 



P/b<l. 



(3.21) 



The first inequality has to be always fulfilled whereas the 
second one has to be satisfied only if the transformation 
is unital, i.e. when f [^/] = hi- Our guess must satisfy 
the condition y/m^ + n? + k"^ < 1 because £[gz] must be 
a density operator. 

In order to proceed further we make an assumption 
that TO = n = 0. By applying the transformation £ ®X 
on the maximally entangled state P+ we have to obtain 
a positive operator. This guarantees the complete posi- 
tivity oi £. In our case, the two smallest eigenvalues of 
the system of two qubits are 



A*± = - ( 1 ± fc 



^ ±2a^cos{e) 



From here we derive two conditions 



l + k> \ oP- 



2q:— cos(0) ; 



l-fc> Wa2 



^ -2a|cos(6'). 



(3.22) 
(3.23) 



If we sum both left and right hand sides of Eqs. (3.22) 
and (3.23), we eliminate the unknown variable k and we 
obtain the condition on the input parameters a, 6, /3 and 
0, namely 



b) - 1 - a2 cos(6')2 ■ 



(3.24) 



The right hand side of Eq. (3.24) is always smaller than or 
equal to unity (equality holds only when the orthogonal- 
ity of the input states is preserved). This is in accordance 
with the requirement (3 < b. Therefore we can conclude 
that the inequality (3.24) represents a stronger condition 
(than /3 < b) that has to be satisfied in order to preserve 
the unitality. 

Because the (pure) state g^ is transformed into the 
state g'z = ■^{I + kS'^), the value of k has to be smaller 
than unity, i.e. k < 1. Under the given condition (3.24) 
the value of k can be chosen to be 



I -a' 



1 — a^ cos^ 



cos 9 < 1 . 



(3.25) 



This value of k fulfills both equations (3.22) and (3.23). 
Consequently, under the condition (3.24) the map £ is 
unital and completely positive. Now our task is to find 
the minimal value of A:. If we try to vary also the pa- 
rameters m, n, then in all numerical tests the distance 
£)(£ [i/],f [^3]) has been found to be larger than the one 
presented above. Consequently, it is preferable to "shift" 
the third state ^3 into a state ^3 which lies on a line 
perpendicular to the plane given by 2^1 ?i 1^2- 

Let us now consider the situation when the total mix- 
ture has to be shifted from the center of the Bloch sphere, 
i.e. our guess is such that £[^I] 7^ ■^I- As we have 
already mentioned above, in this case we have to opti- 
mize over six free parameters. In the previous subsection 
(channel reconstruction based on a measurement of a sin- 
gle test state) the total mixture has been moved along the 
line specified by the center of the Bloch sphere and the 
given state g' . The direct generalization of this feature 
leads us to the following observation: The state £[^1] 
lies in the plane determined by points ^I,£,i,£,'2- In what 
follows we will use the property described at the end of 
the previous paragraph: The third state is mapped into 
a state that belongs to the line perpendicular to the men- 
tioned plane. In other words we set m = n = z = and 
our guess takes the form 



10 

X a^x {(3sm9 — x)/b 

y ~y {(3cos0~y)/b 

fc 



(3.26) 



£ = 



This apporach reduces our task to a three-parametric 
problem which we can solve numerically. We have per- 
formed this task and after finding the solution we were 
searching in the whole six-parametric space for transfor- 
mations that might be better estimates of the map £. 
However no such transformations have been found and 
therefore we conjecture that our estimation is the opti- 
mal one. Unfortunately, we are not able express explicitly 
our guess for general values of a, /3, 6, 9. 




FIG. 4. This figure describes a specific situation of two 
known states £)i,2. Solid lines represent the input "cut" of 
the Bloch sphere and dashed lines correspond to the output 
"cut" of the Bloch sphere. Operators ^[^2 are images of the 
states ^1,2, respectively. Points outside the region of the Bloch 
sphere are associated with negative operators with unit trace. 
The transformation ^j -^ (,'j represents our knowledge about 
the transformation £, i.e. we know that a solid line (given by 
^j) is mapped into a dashed line (given by ^j). In our recon- 
struction scheme the initial states are expressed in Sj basis, 
whereas the output states in the S'j basis. 



1. Examples 



1. Identity. 

Let us assume that the information available about the 
action of the map is of the form Qz ^' Qz = \{I + Sx) 
and ^2 = ^{I + bSy) -^ ^{I + bSy). In this case our guess 
has the form 



/ 1 

10 

10 

Vo fc 



(3.27) 



It is easy to check that the only possibility is to set fc = 1. 
Consequently, our reconstruction is perfect and the iden- 
tity map is uniquely specified. 

The situation is different if we know how the total mix- 
ture is transformed, i.e. f [5/] = ^I- Then our recon- 
struction results in 



10 

10 







(3.28) 



This transformation represents a contraction into the 
line. 

2. Unitary transformation. 

The situation is the same as before. Essentially, we 
have an information about the transformations Qx — *■ 
q'. = W + SL) and 6 = i(/ + bSy) -> i(/ + bS'y). 
Guessing that i/ -^ i/ we obtain Sx — > S'^,Sy -^ S'y. 
Having Sx , Sy then the operator Sz is uniquely identified 
(similarly S'^). The reconstructed map £ has the same 
form as before except it is not defined in a fixed basis. 
This means that our guess is perfect for unitary trans- 
formation. Nevertheless, we have to stress, that the map 
under consideration need not be unitary - it can be a 
non-unital map 

The discussion of the case, when we know that \l ^' 
^I and Qx — *■ q'x is similar to the case of identity. The 
reconstructed map £ has the form (3.28) and again it is 
not expressed in a fixed operator basis. 

3. Contraction to a pure state. 

One can efficiently reconstruct quantum channels of 
this type. The reason is that only when Qx = ^{I+Sx) -^ 
IV') ('01 one is not able to identify such channel perfectly 
with just a single test state. However, if we have knowl- 
edge about two transformations, our reconstruction of 
this channel is complete. 

4. Specific example. 

Let us consider now that we know how two states 

ei-^(/ + 0.6a,); 

g2 = ^(I + 0.4<Jx + Q.l<Jy + 0.8<Jz), 

are transformed under the action of the map (3.11) Ac- 
cording to our strategy we firstly construct two perpen- 
dicular (potentially negative) operators ^1 , ^2 as linear 
combinations of the states gi, g2. One of these operators 
(e.g., <^i) is a pure state. We have two options and one 
of them gives us the pair 

Ci = -(/ + 0.76095(7:^ - 0.080475CTy - 0.6438cr^) ; 

(,2 = ^(1 + 0.464776cr^ + 0.0676122crj, + 0.540897ct;,) . 



We see that the operator ^2 is negative, i.e. it does not 
correspond to any quantum state and Ues outside the 
Bloch sphere. The output states can be expressed as 



From above it follows that the map is not unital and 
therefore, we have to search numerically for the solution. 
As a result we find 



^[ = hl + 0.595857cr^ + 0.34533(7^ - 0.217283cr^) ; 
^^ = -(/ + 0.6402284cra; - 0.0693141^2^ + 0.1825530-^ 



^ad ^ 



1 





0.101568 


0.620589 0.472761 


0.0600669 


-0.0600669 0.620094 


1 


0.457 



and one can easily evaluate the required parameters 

&= 11611=0.71635; 
a= \\(\\\ := 0.722157; 
/3= I |e^2ll =0.669398; 



= ^'^^^"^(^i^'-^^'^ 



== 0.717699 rad. 



Using the two unitary transformations X, Y we obtain 
the guess (the reconstruction of the map £) in the fixed 
basis 



£2^ 



1 











0.1168 


0.83866 


-0.03137 


0.25083 


0.01523 


0.1746 


0.26946 


-0.34023 


0.00696 


-0.0116 


0.36878 


0.28862 



(3.29) 



where we have used the notation ^j = ^{1 + ^j ■ <f)- The 
above formulae follow from the expressions for Ci: '^2^1: '?2 
in Eqs. (3.16) and (3.17), i.e. 



^2 = 1(1 + bSy); 

^2 = 1(1 + (3 sines', + p cos es'y). 



These relations help us to construct the rotation matrices 
X, Y that transform the basis aj ^^ Sj into aj ^^ 5' , and 
vice verse. Our reconstruction results in the matrix writ- 
ten in the adaptable basis, which transforms the matrices 
written in the S'-basis into the 5"-basis. The expression 
of this map in the cr-basis is obtained via the relation 
S2 = Y-'^£'"^X 




FIG. 5. This figure illustrates the transformation of the 
original Bloch sphere under the action of the reconstructed 
map £2 given by Eq. (3.29. The Bloch is transformed into 
an ellipsoid. We can compare this ellipsoid with the ellipsoid 
obtained under the action of the original map presented in 
Fig.(l). 



D. Three input states 



Given three input states we know how a whole plane is 
transformed. Again we have two possibilities: either this 
plane contains the total mixture, or not. If yes, then we 
have the knowledge about the following transformations 

Co = ^/^e^ =l{I + xS', + yS'y); 

ii = lii + s.)^ei = \{i + aS',); 

6 = lil + Sy) ^ ^2 = 1(1 + l3sin0S', +pcoseS'y) . 

Like before, we assume that our experimental data are 
given as Qj -^ g', (j = 1,2,3) and we express the available 
information in a form more suitable for our purposes via 
the states S,j. 

Our aim is to find m, n, k such that the distance 
_D(^Q,^3) is minimal while 



Ca = In 



1 



S.) - C3 = ^{i 



iS' 



nS'^ 



kS'J. 



Of course, the estimated transformation ^ -^ ^3 has to 
fulfill the conditions of the complete positivity. The situ- 
ation is similar as the one discussed in the previous sub- 
section. The only difference is that now we are sure how 
the total mixture is transformed, i.e. whether it stays in 
the center of the Bloch sphere, or not. 

In the second case, when the plane does not contain 
the total mixture, we are forced to make an estimation 
about its position in the state space. It is impossible to 
solve the problem in the most general case for all possible 
parameters. However, in particular cases (having numer- 
ical values of given parameters), we can easily estimate 
the transformation £ in accordance with our strategy. 
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1. Examples 



IV. DISCUSSION AND CONCLUSIONS 



1. Identity. 

In all cases the identity will be identified perfectly. If 
we do not know how the identity is transformed, then the 
guess ^/ -^ \l is reasonable. In fact, we have no other 
choice. Similarly, if we know that -^^ ~* ^^> then our 
only possibility is to set Cs ^ ?3 = Ca- 

We know that the map is unital, i.e. ^7 — *■ ^/ and for 
unital maps with Ai = A2 = 1 the necessary conditions 
are very strict 



1 + fc > 2, 



(3.30) 



and imply that k = \. 

2. Unitary transformation. 

In this case we are dealing essentially with the same 
situation as in the case of the identity. Knowing the 
transformation of three states our guess of the unitary 
transformation is perfect. 

3. Contraction to a pure state. 

As discussed in the previous section this map can be 
completely determined based on the knowledge of how 
two linearly independent states are transformed. 
4-. Specific example. 

Let us assume that transformations of the following 
three states are given (known from the measurement) 

ei = ^(/ + 0.6(7,); 

Q2 = ^{I + 0.4cr, + O.lo-;, + O.Scr^) ; 

Qs = \{I + 0.4a, + 0.3<jy + 0.6a,) . 

Now we are left with only three free parameters and we 
have to search for the reconstruction numerically. The 
result of our numerics is 



/I 

0.294686 0.54219 -0.02396 0.176042 

0.2 -0.3 

V 0.0562144 -0.0936907 0.27918 0.27918 



In this paper we have presented strategies how to re- 
construct (estimate) properties of a quantum channel de- 
scribed by the map £. In a particular case of a qubit 
channel a complete reconstruction of the map £ can be 
performed via complete tomography of four output states 
£[gj] that form at the input of the channel a set of four 
linearly independent states Qj {j = 1,2,3,4). We have 
studied the situation when less than four linearly inde- 
pendent states are transmitted via the channel and mea- 
sured at the output. We have presented strategies how 
to reconstruct the channel when just one, two and three 
states are transmitted via the channel. We have shown 
that unitary transformations (channels) can be uniquely 
reconstructed (determined) based on the information of 
how three properly chosen input states are changed under 
the action of the channel. We conclude the paper with 
three remarks. Firstly we will comment on the channel 
capacity associated with reconstructed (estimated) maps. 
Secondly, we will address the problem of optimality of the 
channel estimation. 



A. The channel capacity. 



Let us study how the channel capacity defined by the 
expression (for details see, e.g. Ref. [6]) 



(3, 




FIG. 6. The picture of the Bloch sphere that is transformed 
by the action of the reconstructed map £3 given by Eq. (3.31). 
Compare with the transformation of the Bloch sphere by the 
map £ in Fig. 1. 



C = max [S{g) 



Hpj'^(^j)] 



(4.1) 



depends on reconstruction strategies. 

Prior any measurement on the channel is performed, 
the conservative assumption about the channel capacity 
should be that it is zero, i.e. the channel is "useless" 
^y the transmission of classical information. One might 
expect that any measurement performed on the channel 
(i.e. via sending a specifically prepared state through 
the channel and measuring it at the output) should re- 
sult in the estimation of the channel that has capacity 
closer to the capacity of the actual channel. However, in 
general this is not always the case. To illustrate this let 
us assume the following problem: Consider that we are 
trying to guess the capacity of the contraction to a pure 
state P^ . If our experiment says that P^ — > P^ , then the 
guessed transformation preserves the total mixture, and 
transforms the whole state space (the Bloch sphere) into 
a line instead into the point. Therefore, the capacity of 
the first estimation is non- vanishing, whereas the actual 
channel has the zero capacity. 
Unital channels 

Let us consider an estimation of a unital channel. The 
the capacity of a unital channel is given by the expression 
[10] 
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C{£) = 1 - H(p) , 



(4.2) 



whereff(A*)--i(l-M)logi(l-Ai)-i(l + Ai)logi(l + 
/i), /i = max{|Ai|, IA2I, IA3I} and A^s are singular values 
of the matrix E. The value of ^ is related to the max- 
imal distance between the total mixture and some state 
g' — £[g\. Because in our estimation strategy we are 
searching for maps that minimize such distance, we al- 
ways reconstruct a channel which has lower (or equal) 
capacity as the actual channel. As a result we obtain 
that for unital maps our estimation strategy is in accor- 
dance with the channel capacity approach: Better our 
estimation is closer is the estimated channel capacity to 
the capacity of the actual channel. 



B. Fidelity of channel estimation 

In order to quantify the fidelity of the channel estima- 
tion we have to introduce a corresponding measure. The 
average fidelity between a map £ and some reference map 
£r can be quantified as an integral 



J^= f d£[d{£,£r)] 



(4.3) 



over all possible maps £. Unfortunately, we do not know 
how to specify a proper integration measure on the space 
of CP maps £. For this reason it is much easier to con- 
sider an average distance between two CP maps £ and Q 
that is defined as 



d{£,g)=j [ 

Js{n) 



[D{£[e\,g[Q\)]dQ- 



(4.4) 



where the average is performed over whole state space of 
the system on which the maps do act. 

A good reconstruction scheme has a property that the 
update of our information cannot debase our estimation. 
In particular, it means that £1 is better estimate as £o: £2 
is better than £1, etc, where £n represents our guess with 
n known state transformations (n = 0,1,2,3,4). Using 
the average distance (4.4) this property can be formally 
expressed via a sequence of inequalities 

d{£n,£) > d{£i,£) > d{£2,£) > d{£s,£) > d{£i,£) , (4.5) 

where £ is the actual completely positive map that is 
estimated and £j are corresponding estimates. 

In our case £0 is a contraction of the whole state space 
into the total mixture. We can evaluate explicitly the 
average distance (4.4) in this case 



di£o,£) 



dg[Tr\ll - £[g]\] 



£is(n)] 



df\r\ , 



where we have used the property D{^I, g) — \r\. Con- 
sequently, d{£o,£) corresponds to a mean distance of 
£[S{H)] to the center of the Bloch sphere. 



Due to the triangle inequality 

D{£n[g],£[g]) + D{£n[g],£i[g]) > D{£i[g],£[g]) 

wc find that 

d{£a,£)+d{£i,£o)>d{£i,£). 

The fact that image of the Bloch sphere under the action 
of the map £1 is a line (a set of measure zero) implies 
that 



d{£o,£i) 



£i[S{H)] 



dr\r\ =0 



and 



d{£o,£) + d{£i,£o) = d{£o,£) > d{£i,£) . 



(4.6) 



From here it follows that the relation d{£o,£) > d{£i,£) 
holds. Consequently the estimation £1 is better than the 
estimation £q. This relation is in accordance with our 
intuition: A guess based on some data must be always 
better than a random guess. 

Following the above argument we conjecture that the 
whole hierarchy of inequalities (4.5) holds. This would 
mean that larger the set of test states better our estima- 
tion is. 

Let us consider the specific example of the map (3.11) 
that has been studied throughout Sec. III. This is a non- 
unital map for which we have presented its estimations 
in various situations, i.e. in cases when one, two or three 
test states have been sent via the channel. Using the 
corresponding estimates £j we can evaluate the average 
distances d{£j,£) for which we find 



(i(£o,f)« 0.845 
(i(fi,£)« 0.758 
d{£2,£)^QMQ 
d{£3,£)KQ.M^. 



We can conclude that the hierarchy of inequalities (4.5) 
for this specific non-unital map is preserved. Once this 
hierarchy is proved to be valid, one can ask a question 
how to chose the set of test states gj so that the sequence 
of distances (4.5) converges to zero (the perfect estima- 
tion) most rapidly. 



C. Optimal test states 

One of our aims was to investigate which states are ef- 
ficient for the process reconstruction. It turns out that it 
is reasonable to start with the total mixture as the first 
test state. Starting with a pure state, the first guess is 
always unital. Using a general mixed state the shift of 
the center of the Bloch sphere is only estimated, but if 
gi ~ ^/, then the question of unitality is solved without 
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any doubts. Therefore, we suggest to use the total mix- 
ture as the first test state. In this case, the channel ca- 
pacity always vanishes after the first step of estimation. 
To improve the reconstruction one has to send via the 
channel more states. The best strategy is to complement 
the total mixture with pure states which arc mutually 
orthogonal in the sense of Bloch sphere representation. 
This optimization of the reconstruction via the choice of 
test states is still an open question. 
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